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J.— s , Both theoretical and experimental results for the dynamics of photoexcited electrons at surfaces of 

("■^ ' Cu and the ferromagnetic transition metals Fe, Co, and Ni are presented. A model for the dynamics 

fvj I of excited electrons is developed, which is based on the Boltzmann equation and includes effects of 

photoexcitation, electron-electron scattering, secondary electrons (cascade and Auger electrons) , and 

transport of excited carriers out of the detection region. From this we determine the time-resolved 
[j' , two-photon photoemission (TR-2PPE). Thus a direct comparison of calculated relaxation times with 

experimental results by means of TR-2PPE becomes possible. The comparison indicates that the 
O^ ' magnitudes of the spin-averaged relaxation time r and of the ratio t^/ti of majority and minority 

relaxation times for the different ferromagnetic transition metals result not only from density-of- 
.j ■ states effects, but also from different Coulomb matrix elements M. Taking Mfc > Mcu > A/ni = 

K^ I Mco we get reasonable agreement with experiments. 
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Q^ ■ I. INTRODUCTION 

on' 

The dynamics of excited electrons at metal surfaces has been studied intensively over the last few years. It is of 

high interest to understand the dynamics of non-equilibrium electrons in different metals on a femtosecond timescale 

and its influence on chemical reactions at surfaces, for example. From the relaxation of hot electrons in ferromagnets 

1-^ , one may also learn about the decay of transient magnetization and about spin-selective transport and tunneling. 

rj ■ Short laser pulses of 15-50 fs duration and the pump-probe technique have made possible the study of electronic 

O . dynamics on ultrashort timescales comparable to typical electron-electron interaction timescales, which in metals are 

O ' of the order of 5-50 fs for excitation energies of about 1-2 eV. 

* * 

h>. The aim of 2PPE experiments is to study the relaxation of single excited electrons. A relaxation time is extracted 

■ ""j from the width of the 2PPE signal as a function of the delay time between pump and probe pulses-S In the earlier 

rN [ experimentsJ^'D the relaxation time was interpreted as the lifetime of a single excitecLelectron due to the Coulomb 

^ ■ interaction. It was compared to the theoretical result from Fermi- liquid theory (FLT).eI The order |e£.njiagnitude and 

■ ■ ■ ' the energy dependence of the FLT lifetime were in good agreement with the experimental results. EIi3'QEI However, it 

was realized that additional physical effects suek as the transport of excited electrons out of the detection region and 

the secondary-electron cascade come into play.Elcl'El It was noted that in 2PPE one generally observes the relaxation of 

a nascent photoexcited electron population and not only the lifetime of a single excited electron. For low excitation 

energies, electronic-lifetimes are longer than a few tens of fs, and ballistic transport leads to a removal of electrons from 

the probed region.Q This is indistinguishable from a stronger electronic decay. Furthermore, under certain conditions 



the measured relaxation time shows a surprising non- monotonous feature which depegids. on the photon energy of the 

m It^was^ 
holes left behind in the excitation can influence the observed relaxation time of hot electrons .EJ^EJ One explanation was 



exciting laser, which cannot be explained by a single-electron lifetime and transport.B tj It pwas. pointed out that the 



that secopdary electrons generated by the filling of holes (Auger electrons) are responsible for the non-monotonous 
behavior.E3 However, the contribution of Auger electrons to the relaxation time has raised some controversy. Petek 
et al. have argued that Auger electrons do not make a significant contributi|OLi|i_tp the observed 2PPE signal and to 
the relaxation time at high intermediate-state energies E — Ep > 1.5 — 2 eV.E3'E3 

The Sd transition metals have not been as intensely studied as the noble metals. However, they offer several 
interesting features which make it worth to study them in detailj-They offer the opportunity to study spin-dependent 
interactions if the spin of the emitted electrons is measured. Ilil'tj Furthermore, the closeness of the d bands to the 



Fermi energy makes it possible to study electronic interactions not only for free-electron-like states, but also for the 
more localized d-electron states. 

On the theoretical side, the effect of the density of states (DOS) on the lifetime [and the influence of secondary elec- 
tronsjoii photoemission from transition metals have been addressed by Pcnn et alx3 Using a similar approach, Zarate 
et aZJ13 have calculated low-energy-electron lifetimes in noble metals and ferromagnetic Co. First-principles liff?time 
calculations have been performed for image-potential statesE^ and for bulk states in alkali and noble metals .E3G3 The 
lifetime is obtained from the inverse of the imaginary part of the self-energy. As in FLT, the lifetime calculated in 
the above works is a single-electron lifetime. Due to the additional effects of secondary electrons and transport in 
2PPE experiments, it is difficult to compare these theoretical results with the relaxation times measured in 2PPE. 
For the 3d transition metals Fe, Co, and Ni, which show important contributions from the more localized d states in 
the vicinity of the Fermi energy, first-principles lifetime calculations in the range of a few eV above the Fermi energy 
have not, to our knowledge, been reported in the literature so far. 

In this paper, we present both theoretical and experimental results for the electron dynamics as observed in 2PPE 
for Cu and ferromagnetic Fe, Co, and Ni. Systematic trends among the transition metals are discussed. A theoretical 
model for the 2PPE process is presented which is based on the time evolution of the distribution function. The 
latter is calculated with the Boltzmann equation including effects of photoexcitation, electron-electron scattering and 
transport. Electron-electron scattering rates are calculated from a golden-rule expression using realistic DOS and 
constat Coulomb matrix elements. The approach for the calculation of the scattering rates is as outlined by Penn 
et alJt3 We extend this approach to include not only the relaxation of excited electrons, but also the generation of 
secondary electrons. Rather than performing a first-principles calculation of the lifetime of single excited electrons, 
we lay emphasis on using a model which yields material-specific single-electron lifetimes for transition metals and 
enables us to calculate the relaxation time of the distribution including effects of secondary electrons and transport. 
This allows a direct comparison of calculated relaxation times with experimental results. 

The structure of the paper is as follows. In Sec. O, we describe the model for the dynamics of excited electrons from 
which 2PPE is calculated. Numerical results for the relaxation of the distribution of excited electrons are presented in 
Sec. III. In Sec. [V, the experiments are described and their results are given. In Sec. 0, experimental and theoretical 
results are compared and discussed. Conclusions and outlook are given in Sec. ^. 

II. THEORY 

The process of two-photon photoemission is illustrated in Fig. [l|. The intensity /^ppe jg obtained by multiplying the. 
distribution function in the intermediate state f{E,a,z,t) with a factor e"^'^" for transmission into the vacuumEll 
and with the power of the laser pulse P{t) and integrating over time t and coordinate z perpendicular to the surface: 

/OO /"OO 

dt P{t) dze-''^^-fiE,cr,z,t) . (1) 

-OO Jo 

Energy and spin of the intermediate state are denoted by E and a. The photon frequency is given by v. For the 
transmission factor, we_use the spin-averaged values of the attenuation length A measured in overlayer experiments 
for different elements.tj The above expression for the photoemission intensity based on the distribution function is 
suited for the description of the population dynamics. Our aim here is to describe incoherent electronic processes 
like the decay of excited electrons and the generation of secondary electrons due to electron-electron scattering. Also, 
the experiments with which we wish to compare our calculations are phase-averaged measurements of the decay of 
the population of excited electrons. The expression is further justified by the fact that for bulk states in metals, 
one expects rapid loss of coherence within a few fs. However, clearly, if one is mainly interested in coherent effects 
like the decay of the optically induced polarization, the treaimept of the dynamics and the photoemission process 
should be based on both occupation function and polarization.aO'Eil Recent interferometric measurements have shown 
relatively long decoherence times in Cu of T^ = 5 — 10 fs for holes at the topLXjf the d bands and electrons at about 
E - Ep = 1 eV and up to T^"^ = 35 fs for electrons at about E - Ep = ^ eV.B 

As shown in Fig. |l|, 2PPE involves three electronic states, in which electron-electron scattering, electronic transport 
and emission into the vacuum take place and determine the observed photoemission signal. After optical excitation, 
the holes left behind in the initial state relax and get filled via Coulomb scattering by electrons from occupied levels 
closer to the Fermi energy. Energy conservation requires that at the same time other electrons from below the Fermi 
energy are excited to unoccupied levels above the Fermi energy (secondary electrons). The holes are also filled via 
transport processes by electrons from the bulk. The optically excited (primary) electrons are scattered out of the 
intermediate state by scattering with electrons in the Fermi sea. On the other hand, secondary electrons are scattered 
into the intermediate state, which leads to the refilling of this state. The intermediate state can be refilled by: i) an 



optically excited (hot) electron after an electron-electron scattering process; ii) a cold electron from below the Fermi 
energy after scattering with a hot electron; iii) an Auger electron (an electron excited from below the Fermi energy 
after a hole is filled by a cold electron) . The latter process leads to a dependence of the observed lifetime on the rate 
of filling of holes (the inverse hole lifetime) . The transport of excited electrons into the bulk leads to the removal of 
electrons from the intermediate state. Thirdly, the final state is above the vacuum energy and describes a free electron 
which can escape from the solid. Only electrons within a mean free path of the surface absorbing a second photon 
are emitted into the vacuum. 

We use now the Boltzmann equation taking into account the above processes to describe the time evolution of the 
electronic distribution function. The electronic states are characterized by energy E, spin a, band index a = sp,d 
and coordinate z perpendicular to the surface. The Boltzmann equation reads 
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and contains the rates of change of the occupation due to optical excitation, electron-electron scattering and electronic 
transport. Note, the relaxation time r of the intermediate state which is compared with the experimental one is 
determined from the decay of the occupation. The details of the procedure are described in App. |^ 

For the ferromagnetic metals, we calculate the relaxation time for excited spin-up and spin-down electrons, T| 
and Tj. Spin-up and spin-down electrons will be denoted as majority and minority electrons in the following. The 
spin-averaged relaxation time r is defined as 1/r — 1/2 (1/t| + l/r|) and the relaxation time ratio as i? = t^ /ti- 

The optical transition rate between two electronic states due to the interaction with the laser field with photon 
energy hv is given by 
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Here, p{Ea, E' (3, v) is an average over electron momenta of the optical transition matrix elements describing the 
transition between an initial occupied state in band a at energy E and a final unoccupied state in band (3 at energy 
E' . The DOS in the final state is denoted by p{E' j3a). We use energy-independent optical transition matrix elements. 
Thus, the strength of the optical excitation is proportional to the initial and final DOS. Optical excitation takes place 
within the optical penetration depth of the surface and has the time dependence of the laser field. We assume weak 
optical excitations, since the energy deposited in the system is not enough to significantly disturb the temperature or 
magnetization. The fraction of excited electrons per atom is about 10~^ (see Sec. IV). 

The transition rates due to electron-electron scattering are calculated using Fermi's golden rule in the random-k ap- 
proximation, since the strong electron-electron interaction in noble and transition metals leads to a fast redistribution 
of electronic momenta so that the information about the initial optical excitation process in k space is quickly lost. 
This then justifies the random-k approximation for the calculation of electronic dynamics. We extend the treatment 
by Penn et a/.E3 to a non-equilibrium situation by calculating scattering rates into and out of a level and taking into 
account the non-equilibrium distribution of electrons. The scattering rates are derived in App. H. The expressions 
for the transition rates for scattering out of or into a state with energy E and spin a =|, | are given by: 
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a,E'a)^— de{e,^h,+^,„2\M^^ + e,sh,+^,^\M^^ 
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and 



W{Ea, E'a) = ^ / deee^he+u^.a |Af ^ ^ | 



(7) 



Here, esa = pEafsa is the number of electrons and hEa — PBa(l — fsa) is the number of holes at energy E, with 
spin a. The energies involved in the transition are E^ E' , e and e + w, where uj = E — E' is the energy transferred in 
the transition. The spin is denoted by a and its opposite by a. For details see App. p. 

To include transport effects, we use a spatially dependent distribution function. Due to the fact that the laser spot 
is much larger than the optical penetration depth, we neglect transport in the direction parallel to the surface and keep 
only the coordinate z in the direction perpendicular to the surface. The distribution function f{E, z) (suppressing 
spin and band indices for simplicity) describes the occupation of an electronic state \E) at the coordinate z . Electrons 
in state \E) are moving in different directions with a certain velocity distribution. In order to describe transport, it is 
necessary to use the velocity in z direction as an additional argument in the distribution function. It is now written as 
f{E, fjfiz)- Using Liouville's theorem, one can express the change in the number of electrons moving with a velocity 
Vr asnj 



df{E,z,v^ 



at 



-v,V, ,fiE,z,v,). (8) 



In the calculation of transport, we take into account the effect of inelastic electron-electron collisions on the distri- 
bution f{E, z, Vz)- In order to do this, we describe how both the electron-electron scattering rates from Eqs. (0-p|) and 
the transport term from Eq. (H) are used in the Boltzmann equation for f{E, z, Vz)- We make the assumption that the 
velocity of the electrons after scattering is randomly distributed. The rate of change oi f{E, z, Vz) due to electrons scat- 
tering into and out of this state is then given by at'^" 
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where the rates on the right hand side are calculated with Eqs. (H-|) using the velocity-averaged distribution 

f{E,z) = (I/A^d) X)i="i fi^T^T^z)- We use Ny discrete velocity intervals of width Avz = {2vp)/Ny), neglect- 
ing the weak energy dependence of the velocity in the range of a few eV around Ep. These terms are used in the 
Boltzmann equation for f{E,z,Vz) together with the terms for transport and optical excitation. The photoexcited 

optical 



electrons have a random distribution of velocities: ot'^' 



"P"'^''' _ df{E,z) 
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The transport effect is caused by the gradient in the particle density created by the photoexcitation within the 
optical penetration depth. We take into account the fact that sp and d electrons have different velocities due to their 
different degree of localization. The velocity of an electron in band a with wave vector k is given by Va{k) — i — ^ ,, ' 
Thus, nearly free electrons in sp-like bands have higher velocities (or smaller effective mass) than more localized 
electrons in flat, d-like bands. For s electrons, we take the Fermi velocity as the maximal velocity. For d electrons, 
we note that the velocity in is roughly proportional to the band_aadth and we therefore set Vd/vs — Wd/Ws- For all 
elements considered, we use vp — 18 A/fs and Wd/Ws — 0.4.lZIH Thus, we distinguish between different elements 
purely by the relative contribution of sp and d electrons. 

We have not included electron-phonon scattering into Eq. (0) because the time scale for energy transfer between 
electrons and phonons is on the order of ps and thus longer than the time scales considered in this work.E3 However, 
electron-phonon scattjesing provides an additional mechanism for momentum transfer and can thus reduce the efficiency 
of ballistic transport H In our calculation, the neglect of electron-phonon scattering might lead to an overestimate 
of the transport effect. The inclusion of the momentum redistribution bi2-.electron-phonon scattering is beyond the 
scope of the present work, but will be the subject of a future publicationHJ 

In Eqs. (y^|5|), the transition rates are determined by the available phase space for a transition weighted by the 
square of the transition matrix element. Scattering out of an excited level and into an excited level are treated on 

the same footing. The two processes occur simultaneously because of energy conservation. The rate L*^ 

for scattering into a level contains the effects leading to the refilling of the intermediate state. By calculating the 
scattering rates in a consistent manner for states above and below the Fermi energy, we keep track of the creation 
and relaxation of electrons as well as holes. 

Note, from the Boltzmann equation we recover the Fermi- liquid behavior t{E) oc (-E — Ep)^'^ for the single- 

Ig^ which yields l/T{Ea) = 1/2 J^;^dE'{hE'aW{Ea,E'a) + 

hE'aW{Ea,E'a)} for the inverse lifetime. For simplicity now we take p = P] = Pi and Af^T — ]\,jM — jvf Then the 
inverse lifetime reduces to 



electron lifetime. To see this we write L °' 
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dE'p{E') de2p{e)p{e + uj)\MY , (9) 
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where co = E — E' . One sees the influence of the DOS within the distance {E — Ep) of the Fermi energy and of the 
Coulomb matrix element. For energies very close to Ep or for constant DOS p one obtains 1/t{E) — ^p^\M\^{E — 



Ep)'^ . Note, that the inverse lifetime is proportional to p^. In addition to the factor p from the relaxation of the initial 
electron at energy E, one also has to take into account the available phase space for the electron-hole pair created 
because of energy conservation, which yields a factor of p^. This may be compared with the FLT lifetime expression 
which is given byEilH t{E) — aQ{rs){E—Ep)~^ , where ao(rs) = 263rr fs eV^ and Ts is the dimensionless parameter 
describing the density of the electron gas. It is given by the relation l/rie ~ 47r(rsao)'^/3, where rig is the electron 

density and uq is the Bohr radius. Thus, the expression for t{E) derived from a*'^ in the appropriate limit 

e— e 

gives the same energy dependence as FLT. 

One can try to understand the relaxation times observed in different metals by means of the FLT expression for 
T by using r^ corresponding to the electron density in the metal. However, the expression is strictly valid only for 
free-electron-like metals and it is not simple to extend it to take into account d electrons. First one can use the 
electron density corresponding to sp electrons only. We take values for integrated sp and d electron densities of states 
in the solid from Ref. p3|. The number of sp electrons per atom for Fe, Co, Ni, and Cu is 1.07, 1.13. 1.03 and 1.10, 
respectively. One obtains roughly r^ ~ 2.6 and qq ~ 25 fs eV'^ for all the metals Fe, Co, Ni, and Cu. This neglects 
d electrons altogether and does not yield any differences between these metals. On the other hand, if one uses the 
total number of sp and d electrons, one obtains values for r^ monotonically decreasing from from r^ — 1.32 in Fe to 
rs — 1.22 in Cu, leading to uq = 130 fs eV^ for Fe and oq = 161 fs eV^ for Cu. It is evident from the magnitude of qq 
that this overestimates the influence of d electrons. Clearly a more refined treatment is necessary which distinguishes 
between sp and d electrons and takes into account the DOS within a few eV of the Fermi energy. The influence of the 
d bands on the lifetime of a state depends on the distance Ed of the d bands to the Fermi energy: for small excitation 
energy AE = E — Ep < Ed, the d bands should have little influence on the lifetime. In Sec. ^, we will discuss the 
influence of d bands in more detail. 

III. NUMERICAL RESULTS OF THE THEORY 

To show clearly the influence of different physical mechanisms, the calculated relaxation times are presented in 
three steps including consecutively more processes in the calculation. 

In the first step, we consider the single-electron lifetime. At this level, we take into account in the Boltzmann 
equation, Eq. (0), only scattering out of a particular level in addition to the optical excitation. Results are labelled 
by (out) in Figs. S and H. After photoexcitation, the time evolution of the distribution function shows an exponential 
decay. The lifetime obtained in this way is a single-electron lifetime and can be compared with FLT or with lifetimes 
obtained from first-principles calculations of the self-energy.EjcJ However, the single-electron lifetime is observed in 
an experiment only if there are no other effects present such as secondary-electron generation or transport. Thus, 
the calculated single-electron lifetime should not be directly compared with experimental results. It serves as a guide 
for comparison with other theoretical results and as a reference for comparison with results when secondary electrons 
and transport are included. 

In the second step, we take into account secondary electrons, while neglecting transport effects. Thus, in the 
Boltzmann equation, we keep the scattering terms for scattering into and out of the levels. Results are labelled {out, 
in). The relaxation time obtained in this way includes effects of the whole distribution. Note, it is not a single-electron 
lifetime, but an effective relaxation time of the distribution. 

In the third step, we also take into account transport effects. The loss of excited electrons due to transport out 
of the surface region will influence the occupation and hence the apparent electronic relaxation time. Results are 
labelled {out, in, transport). 

The DOS for the different—metals used as input in the calculation of the scattering rates in Eqs. ^-^ are taken 
from an FLAPW calculatiorO and are shown in Fig. 0. We distinguish only between d-like and sp-like states in the 
DOS. We take the partial d DOS from the calculation and use the remaining DOS as sp-like DOS. The total DOS are 
very similar to the ones given in Ref. E7l For Cu, however, we shift the d bands to lowep-binding energy by 0.4 eV 
in order to obtain agreement with ARPES results for the binding energy of the d bands.E]l It is known that LAPW 
calculations might yield too small binding energies for the d bands.L3 

The values for the optical penetration depth were obtained from the optical constants in Ref. pfl. We use Apc = 
124 A, Aco — 108 A, Ani — 122 A, and Acu = 149 A for hi^ ~ 3.0 eV. For the range of photon energies considered, 
the energy dependence of the penetration depth can be neglected. 



A. Numerical Results for Cu 

Fig. H illustrates the effect of secondary electrons and transport on the intensity /^ppe (^^ppgi- part) and relaxation 
time (lower part) for optical excitation with photon energy hiy — 3.3 eV. The results for Cu show particularly clearly 
the influence of secondary electrons and transport. We use M = 0.8 eV and |m1^^/M^''-| = 1 for the Coulomb matrix 
elements. This choice will be justified at the end of this section. 

In curve a of Fig. y, we show /^ppe calculated from Eq. (^ without scattering and transport. In this case /^ppe jg 
proportional to the number of optically excited (primary) electrons in the intermediate state. The optical excitation 
in Eq. (^ is proportional to the convolution of initial and final DOS, since we use constant optical transition matrix 
elements. The intensity shows an important contribution from initial states in the d band below E — Ep = 1.3 eV 
with a pronounced peak at 1.1 eV and a small contribution from initial states in the s band above 1.3 eV. 

In curve b of Fig. ^, we show /^ppe keeping only the out-scattering term in Eq. (||), which corresponds to an 
exponential decay of the optically excited distribution. The intensity is reduced compared to the intensity without 
scattering, and the reduction gets stronger towards higher excitation energy due to the shorter lifetime. The relaxation 
time calculated in this way is a single-electron lifetime. For energies below E — Ep = 2 eV, it shows the energy 
dependence t{E) (x {E — Ep)^"^ as in FLT. The FLT lifetime for Cu with cq — 25 fs eV^ is shown by curve e. It is a 
factor of about 2.5 lower than the lifetime calculated using M = 0.8 eV. 

In curve c of Fig. 0, we show the results using the in-scattering term as well as the out-scattering term in Eq. (0). At 
low energy a secondary-electron tail forms and the intensity becomes a superposition of the initial optical excitation and 
the secondary-electron tail. The relaxation time shown by curve c is an effective relaxation time of the distribution 
including secondary-electron effects. It is no longer monotonous and shows a distinct feature in the region of the 
intensity peak. We find a relative minimum at the position of the intensity peak. Further, the relative maximum 
corresponds to the d-band threshold above which very few optically excited d electrons are found. We can understand 
this by studying the contribution of secondary electrons. Comparing the intensities without (curve h) and with (curve 
c) secondary electrons, we find that in the region of the peak, the intensity is only slightly increased by secondary 
electrons, whereas it is strongly increased above the d-band threshold. This shows that in the intensity peak one 
observes mainly relaxation of optically excited electrons, whereas above the threshold, one finds mainly contributions 
of secondary electrons. Comparing the relaxation times, we find that in the region of the intensity peak and also 
at high energies E — Ep > 2.5 eV, the relaxation time with secondary electrons (curve c) comes very close to the 
relaxation time without secondary electrons (curve 6), again showing that one observes mainly relaxation in these 
regions. Most secondary electrons are generated in the process of filling the holes in the d band created by the optical 
excitation (so-called Auger effect). The increase in the relaxation time is due to the fact that secondary electrons are 
generated with a certain delay after the optical excitation corresponding to the d-hole lifetime.Ej'tZl 

In Fig. curve d shows the results obtained using transport in addition to secondary-electron effects. Compared 
to the case without transport (curve c), the intensity is reduced. The transport effect removes excited particles from 
the observation region close to the surface into the bulk and thus reduces the intensity. The relaxation time including 
transport (curve d) has a similar shape as the relaxation time with secondary-electron effects (curve c), but the 
magnitude of the relaxation time is strongly reduced by transport. Interestingly, curve d comes close to curve b in 
the region above the threshold {E — Ep > 1.3 eV). Thus, we find that for Cu in a certain energy range, the effects 
of secondary electrons and traasport on the relaxation time roughly cancel. This has been pointed out in an analysis 
of a 2PPE experiment beforeB3 

Fig. ^ shows theoretical and experimental results for the 2PPE intensity and relaxation time for photon energies 
hv = 3.0 and 3.3 eV. Both the intensity peak and the relaxation time feature shift linearly with photon energy. As 
indicated by the dotted lines, the minimum in the relaxation time corresponds to the peak in the intensity, whereas 
the maximum in the relaxation time corresponds to the d-band threshold. 

We have used a Coulomb matrix element M — 0.8 eV. With this choice we obtain at _E — Ep — 1 eV for the single- 
electron lifetime r = 55 fs and for the relaxation time including secondary electron and transport effects r = 40 fs (for 
photon energy hv — 3.3 eV). i^jhe single-electron lifetime is in good agreement with first-principles calculations in the 
energy range E — Ep < 2 eiZ.E3 It also yields good agreement with experimentally determined relaxation times at low 
energies {E — Ep « 1 eV) .E3 We remark that at higher energies {E — Ep > 2 eV) , we would obtain better agreement 
with first-principles calculations and experiments for a smaller Coulomb matrix element M = 0.6 eV. Thus, the use of 
energy-dependent matrix elements might improve the overall agreement between theory and experiment over a wide 
energy range. 



B. Numerical Results for Fe, Co, and Ni 

First, we discuss numerical results shown in Figs. @a and b and labelled by {out) for the single-electron lifetime in 
Fe, Co, and Ni. We use the same energy-independent Coulomb matrix element M — 0.8 eV and |7\fTT/A/^'''| — 1 for 
the different metals. 

For constant and equal Coulomb matrix element M , the single-electron lifetime is directly related to the DOS shown 
in Fig. g and used as input for the calculation. The influence of the DOS on the lifetime is seen in the scattering- 
rate expression Eq. (0) or in the simplified expression Eq. (0). The scattering rate, the inverse of the lifetime, is 
proportional to a combination of terms which contain products of three factors of the DOS. 

When comparing the single-electron lifetime of Cu shown in Fig. (curve b) with the results in Fig. Ha for the 
transition metals, one can see that Cu has a lifetime much longer than the other metals. This is due to the small total 
DOS close to Ef- In Cu, the d bands are located about 2 eV below Ep and there is only a very small total DOS close 
to Ep (Fig. H). The small d DOS close to Ep is due to hybridization with sp-like states. In Ni with one electron less 
than Cu, the d bands move closer to Ep. Furthermore, the d bands are split into a minority and majority spin band 
and a small portion of the minority d band is unoccupied, extending up to about 0.4 eV above Ep (Fig. ||). Due to 
the pronounced peak at the upper edge of the d-band, the minority DOS close to Ep is extremely large. This leads 
to a large phase space for electron-electron scattering at low energy. Thus at low energy (below E — Ep = 1 eV), 
Ni has the smallest calculated single-electron lifetime among the four elements. Co, with one electron less than Ni, 
has an even larger portion of unoccupied minority DOS, extending up to about 1.2 eV above Ep. Although the total 
number of unoccupied states is higher in Co than in Ni, the minority DOS at Ep is lower in Co. Thus at low energy 
(below 1 eV), Co has less phase space and the calculated lifetime is longer than in Ni. With increasing energy, more 
and more unoccupied states become available in Co, so that above 1 eV, the calculated lifetime in Co becomes shorter 
than the one in Ni. In Fe, again with one electron less compared to Co, the unoccupied minority DOS extends up 
to 2.4 eV above Ep, and even the majority DOS has a small unoccupied fraction. The minority DOS at Ep in Fe is 
lower than in Co and in Ni, so that Fe has the smallest phase space and the longest calculated lifetime at low energy. 

Therefore, at low energy (below E — Ep = 1 eV), our simplified theory with equal AI for the different metals gives 
^Ni < Tco < Tpe- This trend changes for Co and Ni above 1 eV. Then we get rco < tni- At even higher energy (above 
2 eV, not shown in the figure), all of the unoccupied d states in Fe are available for a transition and one calculates 
TFe < Tco < TNi- This relation is also, observed in transmission experiments above the vacuum energy for electrons 
with energies above E — Ep = 5 eV.E3 

In Fig. ^ the ratio of majority to minority single-electron lifetime t^/ti is shown. For Co, this ratio is nearly 
constant with a value t-^ /t^ — 7.5. This is understandable in view of the high and nearly constant ratio of minority 
to majority DOS at low energy. For Ni the ratio decreases from t-^/ti = 9.5 a.t E — Ep = 0.4 eV to r-f/rj^ = 4 at 
E — Ep = 1.4 eV. The decrease is due to the fact that above 0.4 eV there are no more unoccupied minority d states. 
The additional phase space gained by going to higher energy is the same for minority and majority electrons, leading 
to a smaller ratio. In Fe the ratio increases from t^/ti = 0.5 to 1 for excitation energies between E — Ep = 0.4 and 
1.2 eV. Thus, majority electrons have a shorter calculated lifetime than minority electrons at low energy. This results 
from the unoccupied portion of the majority DOS above Ep, which for low energy leads to a larger phase space for 
the relaxation of majority electrons and therefore to a shorter lifetime. 

Secondly, we discuss results shown in Figs. 0c and d obtained when secondary electrons are included in the cal- 
culation. They are labelled by {out, in). Transport effects are still neglected. The inclusion of secondary electrons 
leads to an increase of the relaxation time by a factor of about two as compared to the single-electron lifetime. The 
strongest effect is found at the lowest energies. The increase is stronger for the elements with the shortest calculated 
lifetimes (Ni, Co), so that the differences in the relaxation time including secondary electrons between Ni, Co, and Fe 
are smaller than for the single-electron lifetime. The ratio t^ /ti for Ni and Co is reduced to t^/ti = 4 — 5 for Ni and 
t^/t^ =5 — 6 for Co. For Fe, the ratio is nearly unchanged, t-^/ti — 0.5 — 1. The reduction of the ratio t-^ /ti due to 
secondary electrons is understandable in view of the fact that the inclusion of secondary electrons leads to a coupling 
between majority and minority electron populations via electron-electron scattering. Relaxing minority electrons can 
excite majority electrons and vice versa. For Co and Ni, for example, longer-living majority electrons will continue 
to excite minority electrons after the shorter- living primary minority electrons have relaxed. The apparent minority- 
electron relaxation time becomes longer and the ratio tt:/ti becomes smaller by this process. The trends among the 
calculated relaxation times of the transition metals, particularly the relation tni < tco < tfo at .E — Ep < 1 eV are 
unchanged when secondary electrons are included. 

The spectral shape of the optical excitation, i.e. the distribution of primary electrons, has some influence on the 
calculated relaxation time. For example, electrons excited to a high energy lead to more secondary electrons (due 
to the short lifetime of the primary electrons) and to a distribution extending to higher energy (due to the high 
energy of the primary electrons). Similar arguments apply to the energetic position of the holes created by the optical 



excitation. The results for the transition metals shown in Fig. ^ are obtained using an optical excitation with photon 
energy hv — 3.0 eV and constant optical transition matrix elements in Eq. (g). We have also studied the effect of a 
different optical excitation on the relaxation time. We have modelled a resonant optical excitation where excitations 
take place dominantly between an initial state at the top of the d band and a final state in the sp band. No significant 
difference in the calculated relaxation time for the two different shapes of the excitation is obtained. 

Thirdly, we discuss results of the calculations including transport in addition to secondary-electron effects in the 
Boltzmann equation, Eq. (0). They are labelled by {out, in, transport). Results are shown in Figs, rae and f. The 
inclusion of transport leads only to a very slight reduction of the effective relaxation time in Fe, Co, and Ni. The 
ratio T^ /ti as well is only slightly affected by transport. This is in contrast to the strong reduction of the effective 
relaxation time in Cu (c»pi|pare curves c and d in Fig. y). The explanation is that the typical transport relaxation 
timescale is about 40 fs.alij This is shorter than, or comparable to, the single-electron lifetime in Cu in the energy 
range considered, but considerably longer than the typical lifetimes in Fe, Co, or Ni. Thus, transport is expected to 
have great influence for Cu, but not for the transition metals. 

IV. EXPERIMENT 

A. TR-2PPE technique 

The TR-2PPE pump-probe experiments are carried out in a UHV chamber by monitoring the number of electrons 
at a given kinetic energy as a function of the delay between the pump and probe pulses. We employ the equal- 
pulse correlation technique, i.e. the two pulses are monochromatic and equal in intensity, but cross-polarized. For 
metals, the use of orthogonal linear japlarized light pulses suppresses coherent interference effects (within the limit 
of rapid dephasing) to a large extent u Otherwise they influence the optical transition process and would make the 
reconvolution of the raw data much more difficult. Furthermore, the influence of Cs-induced surface states on the 
lifetime can be suppressed (see below). 

The non-linearity of the two-photon process leads to an increase in the 2PPE yield when the pulses are spatially 
and temporarily superimposed. As long as the two laser pulses temporarily overlap it is obvious that an electron can 
be emitted by absorbing just one photon from each pulse. However, if the pulses are temporarily separated, then 
an excited electron from the first pulse is able to absorb a photon from the second pulse but only as long as the 
inelastic lifetime of the intermediate state exceeds the delay or the normally unoccupied electronic state is refilled by 
a secondary electron. Due to a precise measurement of the time delay between the two pulses (1 fs = path length 
difference of 0.3 /xm), this technique allows us to analyze relaxation times which are considerably shorter than the 
laser pulse duration. 

We use laser pulses at low fluence and peak power to avoid space-charge effects or highly excited electron dis- 
tributions. We emphasize that the count rate is much lower than one electron per pulse. Therefore, we measure 
the relaxation and transport of individual excited electronic states rather than the collective behavior of transiently 
heated non-equilibrium distribution. We have to roughly calculate the fraction of excited electrons. Typically, we 
have a laser fluence of about 0.3 nJ/pulse in each beam resulting in 6 x 10* photons per pulse. For a spot size of 
~ 150 /xm and a penetration depth of the blue light of ~ 150 A, the volume in which the laser light will be absorbed 
is about 3 x 10"^" cm'^. If 10% of the light is absorbed by the metal, then 6 x 10^ photons are absorbed by 7 x 10^"^ 
atoms which results in a fractional excitation of roughly 1 part in 10^. 

B. Experimental set-up 

A schematic overview of the experimental set-up is shown in Fig. a The samples are mounted in a UHV chamber 
with a base pressure in the 10^^^ mbar range. It is equipped with a cylindrical sector electron. energy analyzer (CSA) 
and a spin analyzer, based on spin-polarized low-energy electron diffraction (SPLEED).l3 The earth's magnetic 
field is shielded by /i-metal coverings inside the chamber. Standard surface-physics methods such as Auger-electron 
spectroscopy (AES) and low-energy electron diffraction (LEED) are available to check the cleanliness and the surface 
structure of the samples. The orientation of the samples is 45° with respect to the laser beam and the electrons 
are detected in normal-emission geometry. Remanent magnetization of the ferromagnetic samples is achieved by a 
magnetic field pulse from a coil. The geometric arrangement of the spin analyzer allows the measurement of the spin 
polarization along the horizontal in-plane axis of the sample. To minimize the effects of stray fields and to facilitate 
electron collection, a bias voltage (—4 V for Cu and —15 V for Ni, Co, and Fe) is applied between the sample and the 
CSA. For spin-resolved measurements the electrons are guided into the SPLEED analyzer, which is located on top of 



the CSA. In this spin analyzer the electrons are first accelerated to 104.5 eV-ikinetic energy, as the highest figure of 
merit for this kind of analyzer is known to be at this primary-electron energy.Cj Thereafter, the electrons are scattered 
at a tungsten (001) crystal. From the resulting LEED feature, the (-2,0) and (2,0) electron beams that have a high 
spin asymmetry, are counted in two different channeltrons. The Sherman factor S, a quantity for the spin selectivity 
of an analyzer, is found to increase from 0.2 to 0.25 over the two years of operation; the highest value was reached 
after having the tungsten SPLEED crystal a long period in very low pressure. 

The time-resolved 2PPE experiments are performed with a femtosecond mode-locked Ti:sapphire laser, pumped by 
about 10 W from a cw Ar"*" laser. The system delivers transform-limited and sech^ temporal shaped pulses of up to 9 
nJ/pulse with a duration of 40 fs at a repetition rate of 82 MHz. The hnearly polarized output of the Ti:Sapphire laser 
is frequency-doubled in a 0.2 mm thick Beta Barium Borate (BBO) crystal to produce UV pulses at /ii/ = 3 to 3.4 
eV. The UV beam is sent through a pair of prisms to pre-compensate for pulse broadening due to dispersive elements 
like lenses, beamsplitters and the UHV-chamber window in the optical path. A GVD and intensity-loss matched 
interferometric autocorrelator set-up is used for the pump-probe experiment (see Fig. 0). The pulses are split by 
a beamsplitter to equal intensity (pump and probe pulses), and one path is delayed with respect to the other by a 
computer-controlled delay stage. Both beams are combined co-linearly but cross-polarized by a second beamsplitter 
and are focused at the sample surface. 

For the ferromagnetic samples, we use evaporated films because they can be held magnetized in a single-domain 
state without an applied external field and the stray field is much smaller compared to a bulk ferromagnet. In 
principle, this experiment could also be performed with bulk samples. The ferromagnetic films are evaporated onto a 
Cu(OOl) substrate in a separable chamber. We use a water-cooled evaporator based on electron-beam heating. The 
material to be evaporated (of 99.999% purity) was inside a molybdenum crucible (Co, Fe) or directly evaporated from 
a 1 mm thick wire (Ni). The evaporation rate (around 0.2 nm/min) was checked with a quartz oscillator, which is 
calibrated against atomic-force-microscope thickness measurements. During evaporation the pressure remained in the 
10"^" mbar region. 

The thickness of the ferromagnetic films must fulfill the following requirements. Firstly, it has to be large enough, 
in order to avoid an infiuence of electrons from the Cu substrate. Secondly, it should be possible to remanently 
magnetize the film by a suitable strong field pulse. And thirdly, the axis of the magnetization has to lie in the film 
plane because the geometry of our spin analyzer only allows the measurement of transvcrsally polarized electrons. 

Cobalt films, epitaxially grown on Cu(OOl) surfaces, are formed in a stable fee structure and exhibit in-plaii£ 
magnetization. For thick films, the in-plane magnetization easy axis lies along the (110) direction of the Cu crystal.cll 
In-plane magnetization was detectable starting at a_thickness of around 0.4 nm. This is in agreement with other 
investigations using the magneto-optical Kerr effecto Above a film thickness of 2 nm, the spin polarization did not 
increase any more on further deposition of Co. For our investigations we evaporated 10 nm thick Co films. 

Iron grows in the fee structure on Cu(OOl) during the initial steps of evaporation. The magnetization vector is. 
oriented perpendicular to the film surface. Above a thickness of around 2 nm a bcc (110) steucture starts developing .c2l 
The magnetic easy axis for these thick films is found to lie parallel to the Cu(lOO) axisa We use 20 nm thick iron 
films for our investigation. 

For Ni on Cu(OOl) the niagnetization vector is in-plane for small thicknesses, then it switches to out-oLplane at a 
thickness of around 1.2 nmo Only at larger thicknesses of around 6-7 nm does it turn back to in-plane.c3 We found 
a saturation of the spin polarization for thicknesses above 25 nm. Therefore, we evaporated 40 nm thick Ni films for 
our measurements. 

The clean metal surfaces are first dosed with Cs to lower the surface work function, a well-known technique. This 
enabled lifetime measurements of lower excited states, increasing the useful energy range of the spectra closer to the 
Fermi energy (see Fig. |l|). Cs is evaporated from a thoroughly out-gassed getter source (SAES). The effect on the 
lifetinjie.|tiy dosing a metal surface with small amount of Cs (<0.1 ML) has been thoroughly investigated in the last 
years.cZlc3 Using cross-polarized pulses no di^rences in the lifetime have ever been found between a clean and a 
cesiated metal surface by means of TR-2PPEH3 In addition, we found no differences in spin polarization between the 
clean and the cesiated surfaces in the overlapping energy region between 1.7 and 3.3 eV. 

C. Experimental results for Cu 

In Fig. Hb the extracted relaxation time as a function of the intermediate-state energy for a Cu(lll) surface is 
shown, using a photon energy of 3.0 eV (□) and 3.3 eV (O). The data are reconvoluted from the experimentally 
obtained cross-correlation traces using a rate-equation model for the population of the intermediate state. In the 
case of rapid dephasing, and assuming an exponential depletion of the nascent photoexcited electron population, 
the evolution of the transient population N*{t) of the intermediate state is given by dN*{t)/dt = A{t) — N*{t)/T 



where A{t) is the excitation induced by the first (pump) laser pulse. Ekjtails of the way to extract the relaxation 
time from the measured signal have been given in previous publications .tJu'cHI As expected, the lifetime increases as 
the e xcited state energy decreases, caused by the reduced phase space for electron-electron scattering (see Sees. ||, 



III A , and Fermi- liquid theory). However, at an intermediate-state energy, where the intensity becomes dominated by 
intcrband transition from the d band to the unoccupied sp band (strong rf-band peak in the intensity), the measured 
relaxation time decreases by more than a factor of two before it increases again. By changing the photon energy from 
hh' — 3.0 eV to hh' — 3.3 eV, both the d-band peak in the intensity and the dip in the measured lifetime move with 
the same energy difference AE = Ahv as expected. The striking difference in the values observed in the energy range 
E — Ep = 0.7 — 1.3 eV indicate quite clearly that the relaxation time can depend critically on the used photon energy. 

D. Experimental results for Fe, Co, and Ni 

1. Spin-integrated time-resolved 2PPE measurements 

We used the same equal-pulse correlation technique to extend the investigation of the hot-electron relaxation to 
transition metals Co, Fe, and Ni. Compared with noble metals, in which the d shell is completely filled (see Cu and 
Ag), the d band of transition metals is only partially filled, and the electronic and relaxation properties are dominated 
to a considerable degree by these d electrons. The strong localization of these d electrons results in a narrower band 
and hence in a much higher DOS near the Fermi level as compared with the sp electrons in Cu and Ag. A higher 
density of occupied and unoccupied states near the Fermi level is expected to lead to faster relaxation and hence to a 



shorter inelastic lifetime of excited electronic states as discussed in Sec. IIIB. This prediction is well satisfied by our 
data. Fig. [^ shows a comparison of the extracted relaxation time of silver and the three investigated ferromagnetic 
transition metals cobalt, nickel and iron. The experimental values for these metals are at least a factor of 10 smaller 
than those of Cu and Ag. These small values reduce the energy range in the region of 0.3 eV to 1.3 eV which provides 
a meaningful statement about the relation in the relaxation time between the three investigated transition metals. 
In contrast to the numerical results calculated using the same Coulomb matrix element M for different metals (see 



Sec. [II B), the data indicate a relation rpe < TTsii < tco within this energy range. 



2. Spin- and time-resolved 2PPE measurements 

Adding a spin analyzer to the CSA energy analyzer makes possible the separate but simultaneous measurement of 
both spin states. The electrons of a fixed energy are counted according to their spin in two different channeltrons 
as a function of the time delay between the two pulses, at a given magnetization direction. To compensate for an 
apparatus-induced asymmetry, the magnetization is then reversed and the measurement is taken again. From the 
resulting four datasets the relaxation times T| and r^ for spin-up and spin-down electrons are extracted by using the 
same reconvolution method as discussed above. Each pair of data points presented in the plots of this section is the 
average of eight to ten single relaxation time measurements. 

The spin dependence will have a superimposed effect on the spin-integrated relaxation time. Therefore, a spin 
dependence in the relaxation time can only be resolved if there is already a certain relaxation time detectable with 
spin-integrated measurements. As shown in Fig. R in the energy range above 1.4 eV, we find for all three transition 
metals a relaxation time smaller than our time resolution (<2 fs). On the other hand, at intermediate-state energies 
close to Epj the electrons emitted by IPPE processes start becoming important. They induce a large background 
to the 2PPE signal and make an accurate extraction of the lifetimes difficult. Therefore, spin-resolved measurements 
can only be usefully performed for intermediate-state energies between 0.3 and 1.1 eV. 

In Fig. ra the spin-dependent relaxation time for electrons (upper part) and the ratio of majority to minority 
lifetime (lower part) of Fe, Co, and Ni films are plotted. The error bars in the plot represent the statistical scatter. 
The experimental results of the three examined ferromagnetic materials show two common facts: i) The lifetime for 
majority-spin electrons is always found to be longer than the lifetime for minority-spin electrons and ii) the value for 
ti^/ti was found to lie between 1 and 2. The largest differences between T| and T| are found for Ni and Co, whereas for 
Fe, the difference is slightly reduced. This qualitative behavior of the spin-dependent lifetime can be readily explained 
by the excess of unfilled minority-spin states compared to unfilled majority-spin states. According to this simple 
model, the spin dependence of the scattering rate is larger for the strong ferromagnets Co and Ni than for the weak 
ferromagnet Fe. This is in agreement with our measurements, where only a small spin effect could be detected for Fe. 
In Fe, this model would even predict a reversal of the effect for low energies below 1 eV, i.e. the lifetime for spin-down 
electrons should become longer than the lifetime for spin-up electrons, see Fig. ob. A ratio of majority to minority 
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relaxation time R below 1 is, however, not observed for E — Ep < 1 eV. This result indicates that the simple model, 
considering the different number of empty electronic states as the only decisive factor for a spin-dependent relaxation 
time, is not sufficient for a quantitative interpretation of our experimental data. 

V. DISCUSSION 

First, in Fig. ^ we compare experimental and theoretical results for Cu for photon energies hv = 3.0 and 3.3 eV. 
Experimental and theoretical results show qualitative agreement regarding the main features. The peak in the inten- 
sity and the feature in the relaxation time (relative minimum and maximum) shift linearly with photon energy. The 
minimum in the relaxation time corresponds to the peak in the intensity and the maximum in the relaxation time 
corresponds to the d-band threshold. The ex planat ion that the feature in the calculated relaxation time is due to the 



secondary electrons was given in detail in Sec. |III A| . The good agreement between calculated and experimental reswlte 
is a strong evidence for this explanation. The calculations reproduce the features observed in different experimentsETta 
in a natural way by|-including secondary electrons without invoking further explanations such as excitonic states in- 
volving 3d electronsJa The differences between theoretical and experimental results in Fig. ^lead to several conclusions. 
Both the peak in the intensity and the difference between the minimum and maximum of the relaxation time are 
more pronounced in the experiment than in the calculation. This may be an evidence that the calculation yields too 
many secondary electrons which cause a too strong background and thus a too small structure in the 2PPE intensity 
(see the upper part of Fig. Ha). The relatively small increase from the minimum to the maximum in the calculated 
relaxation time should not be affected much by this, but may rather point to the fact that the lifetime of the d holes 
in the calculation is too small. Note, due to the large d DOS below Ep, there is no symmetry between hole lifetimes 
and lifetimes of excited electrons. In our calculations, hole lifetimes for energies below the d-band threshold are very 
short due to the large DOS. If the hole lifetimes were larger, then fewer secondary electrons would be generated, but 
with a longer delay after the creation of the electron-hole pair by the laser pulse. A longer delay will lead to the 
observation of a longer relaxation time in the excited state when secondary-electron contributions are important. 

The interpretation given here that the non-monotonous-feature in the lifetime iOnCu is due to secondary (Auger) 
electrons has raised some controversy in the literature. Bild'E3E3 While Knoesel et alS3 interpret-their data by contribu- 
tions from Auger electrons at intermediate-state energies above the d-band peak, Petek et alS3 argue that secondary 
electrons make no significant contribution to the signal above E — Ep = 1.5 eV (for photon energy hi/ = 3.1 eV). The 
argument is based on a temperature-dependent delayed rise in the 2PPE signal as a function of time delay between 
the laser pulses, which is observed below and immediately above the d-band peak a,t E — Ep = 0.9 eV, but which 
vanishes above E — Ep = 1.5 eV and in the region of the peak. The temperature-dependent delayed rise is interpreted 
as a contribution from Auger electrons, which agrees with our interpretation of the feature in the lifetime. In Ref. 
O, the fact that the delayed rise vanishes above 1.5 eV is taken as evidence that Auger electrons are absent in this 
energy region. In contrast to this conclusion, our calculations show significant contributions from secondary electrons 
up to about E ~ Ep = 2.5 eV (compare the relaxation times without [curve 6] and with [curve c] secondary electrons 
in Fig. 0). We argue in the following that the absence of a resolvable second peak in the 2PPE signal is no evidence 
for the absence of Auger electrons. Thus the results of Ref. |l3| are not in contradiction with our results. One would 
observe a second peak with a delay given by the hole lifetime at the d-band peak if all Auger electrons were created 
at a fixed rate corresponding to this hole lifetime. However, Auger electrons are also created by the filling of holes 
deeper in the d band with energies up to hv. These deep holes have have shorter lifetimes than the ones at the top 
of the d band. Thus they lead to secondary-electron contributions to the dynamics in the intermediate state with a 
smaller delay time. The fact that holes with different lifetimes contribute to the secondary-electron dynamics makes 
it difficult to observe a resolvable second peak with a fixed delay corresponding to the lifetime at the d-band peak. 
Thus in our view, the measurements reported in Refs. P^,f3^ are not in contradiction with the interpretation of the 
non-monotonous feature in the relaxation time given by us and in Ref. HOl 

In Figs. and M we compare experimental and theoretical results for the spin-averaged relaxation time r and the 
ratio T^^JT^ of majority and minority relaxation time for the ferromagnetic transition metals Fe, Co, and Ni. The 
discrepancies between experimental and theoretical results indicate that both DOS and Coulomb matrix elements 
play a role. Note, theoretical results refer to relaxation times of the distribution including secondary-electron effects. 
Transport effects have been neglected here in view of the fact that they cause only minor changes in the relaxation 
time of the transition metals, see Fig. pi 

First, the results calculated with M = 0.8 eV for all the transition metals are shown by the curves a in Figs. ^ and 
QQ. The difference in the calculated relaxation times for Fe, Co, Ni, and Cu is then only due to the different DOS 
used as input for the calculation. Note, the calculated relaxation time is smaller than the experimental one in Co 
and Ni, while it is larger in Fe. The calculated ratio ti/t^ is larger in Co and Ni than the experimenal one, but it is 
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smaller in Fe. 

Secondly, in curves b, we show results of calculations using again M = 0.8 eV, but the reduced value \M^^ /M^^ | = 
0.5. One expects that the matrix element M'^^ for scattering of parallel spins is smaller than M^^ for antiparallel spins, 
since electrons with parallel spins avoid each other due to the Pauli exclusion principle.c2l The ratio t^ /ti is strongly 
reduced in Co and Ni, while it is increased in Fe, which leads to satisfactory agreement for t^/ti between experimental 
and theoretical results. The spin-averaged relaxation time is not strongly affected by the value of |MTT/MT'''|. 

Thirdly, we take into account different Coulomb matrix elements M for the various metals, while we still use 
|mTT/m1^-''| — 0.5. The results are given by the curves c in Figs. and [iQ. For Co and Ni we use M — 0.4 eV, while 
for Fe we take Af = 1.0 eV. The use of these values for M leads to reasonable agreement between theoretical and 
experimental results for both the spin-averaged relaxation time and the ratio t^/t^. 

Different Coulomb matrix elements in Fe, Co, Ni, and Cu are mainly caused by-tiie influence of d electrons. Note, 
while in isolated atoms. Coulomb matrix elements do not vary much from Cu to Fe,Lll in solids the band character, the 
position of the d band, and the screening of d electrons are expected to change this. The screened Coulomb matrix 
elements for scattering between Bloch states with wave vectors k^ in bands a^ is given by 

M's;! = J d\ d\' C(r)Cc..(rO ^(|^_^,|'^)|^_^,| ^k3a3(r)V^k.o.(rO • (10) 

Here, e is the dielectric function. The Coulomb matrix elements are of course influenced by d electrons, since their 
wave functions are more localized and also since they contribute to the screening of the Coulomb potential. The 
strong localization of d electrons leads to smaller overlap with sp-electron wave functions and therefore to smaller 
transition matrix elements when sp — > d transitions are involved as compared to matrix elements involving sp — > sp 
transitions. Note, the d-electron wave functions get more localized from Fe to Cu. The additional screening of d 
electrons is contained in the dielectric function £(|r — r'|,ti;), where lo is the energy transfcred in the transition. 
Somewhat depending on w, d electrons closer to the Fermi energy contribute mainly to screening. In the static limit 
(w ^- 0), the Lindhard dielectric function for a free electron gas reduces to e{q) = 1-1- k^/q^ such that the screened 
Coulomb interaction in real space takes the form V{r) — ^e~Ji°'^ . In the Thomas-Fermi approximation, the screening 
wave vector is directly related to the DOS at the Fermi level,E3 fco = A'Ke'^p{Ep). In the case of transition metals, the 
expressions are not strictly valid because d electrons are not free-electron-like. Although the quantitative contribution 
of d electrons to screening is not well-known, qualitatively it is clear that a higher DOS near the Fermi level leads to 
stronger screening. This may explain that the screened Coulomb matrix element in Co and Ni with many d electrons 
close to the Fermi energy is smaller than the one in Cu with nearly no d electrons close to the Fermi energy. For 
larger energy transfer w, also electrons further away from the Fermi energy contribute to screening. Then ultimately 
the total number of d electrons influences screening. This could be the reason why Fe, which has fewer d electrons, 
has a larger Coulomb matrix element than Co, Ni, and Cu. 

After completion of our study we became aware of ^•platprl work about the influence of d electrons on the lifetime 
of low-energiL-clectrons in noble and transition metalsO'Oo Choosing [different matrix elements for sp and cLatates, 
Zarate et a/.EII obtain good agreement with experimental results for Co.EJ Campillo et a/.Ej and Schone et al\^ have 
calculated lifetimes in Cu using the density-functional theory for the determination of the electronic structure and 
have found important contributions of d electrons to the lifetime via screening, localization of the wave function and 
DOS effects. 

VI. CONCLUSION 

We have presented experimental and theoretical results for the dynamics of excited electrons in Cu, Fe, Co, and Ni. 

The results for Cu show the influence of secondary electrons and transport effects on the observed relaxation 
time. The non-monotonous behavior in the relaxation time obtained in the calculation is in qualitative agreement 
with experiments. It seems desirable to achieve a better quantitative agreement for this structure in order to draw 
definitive conclusions about the role played by secondary electrons. 

Experimental results for the spin-dependent relaxation time for Fe, Co, and Ni reveal that rpo < tni < tqo and 
that t-\^/ti lies between 1 and 2 for the three metals. 

The comparison of experimental and theoretical results shows that DOS effects alone do not explain the magnitude of 
the relaxation time for the various transition metals observed in experiments. The differences between the calculation 
using the same Coulomb matrix element for Fe, Co, and Ni and experiments reveal that Coulomb-matrix-element 
effects are important. 

As an outlook for further studies, we conclude that more detailed calculations have to include a first-principles 
calculation of the Coulomb interaction matrix elements. Notably one has to take into account the screening by 
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the d electrons, the influence of the localized d-electron wave functions on the Coulomb matrix elements, and the 
energy dependence of the matrix elements. Further studies of the influence of the optical excitation, for example the 
influence of the photon energy and the polarization of the incoming light on the observed dynamics, are needed. Also 
the influence of the hole lifetime on the calculated relaxation times should be investigated. The rate of filling of d-band 
holes influences the time evolution of the distribution and hence also the observed relaxation time for levels above the 
Ferpii energy via the generation of secondary electrons. Hole lifetimes have been observed in recent expftcipients on 
Cu|^ and their effect on two-photon photoemission has also been considered in a recent theoretical work.ll2l 
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APPENDIX A: DETERMINATION OF THE RELAXATION TIME 

In order to determine the effective relaxation time, we fit the occupation function calculated from Eq. (0) with a 

optical 



function / describing exponential decay with a relaxation time r obtained from the equation: 7^ = 7^ 



/ 



dt ~ dt 

The fit to the occupation calculated from Eq. (g) is done by taking r for which / has the maximum at the same 
time. The comparison of the curves in Fig. y_l| shows that this is a valid procedure over a wide range of energies, even 
where secondary electrons dominate (see E — Ep ~ 0.8 or 1.3 eV). The fact that the deviations are very small shows 
that although the calculated curves do not exactly show exponential behavior, they can be fitted well by a curve / 
showing exponential decay for some effective relaxation time. 

APPENDIX B: ELECTRON-ELECTRON SCATTERING RATES 

The scattering rate out of the state with momentum p, band a (designating sp or d-like wave function) and spin a 
is given in first order time-dependent perturbation theory (golden rule) by: 



(y Jpaa 



dt 



out 



— -./pacr-r- 2^ S -^fkyai^ - /p'/3cr)(l ^ fk'Sa) 



p'kk',/37(5 
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In the same manner, one defines the scattering rate g°' 



for scattering into a state. The first and second terms 



describe scattering between electrons of of the same and of opposite spin, respectively. ■— ■ 

Sums over momenta are converted into integrals over energies in the random-k approximation.t^ The conversion to 
k-averaged quantities is done in the following way: 

/oo 
dEpEaafsaa- ■ 
K -°° 

Each k sum leads to a factor of the DOS, psaa- Products of distribution functions and densities of states are defined 
as CEcxa — PEaafEaa and hEaa = PEaa{^ — /bqct) and designate the number of electrons and holes. The random-k 
approximation allows the replacements p ^ E, p' -^ £", k — > e and k' -^ e' . The integral over e' can be performed 
because of the (S- function and allows to replace e' by e -I- oj with uj = E — E' . The k-averaged expressions for the 
scattering rates out of and into the state Eaa are then given by: 

13 
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For the calculation, we will use a simple parametrization of the matrix elements. Fi rst, we do not further distinguish 
between states of sp or d symmetry in the matrix element. Then the form of Eqs. (B2-B5) remains the same when 
the partial DOS is replaced by the total DOS. Secondly, we neglect the energy dependence, assuming it to be weak 
in the range of a few eV from the Fermi energy. In Eq. (B4), we neglect the interference term after expanding the 
modulus square: 
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(B6) 



We have denoted the matrix elements with different energy arguments by A/f and M|''^. In the last step, we 
drop the energy index and use energy- independent matrix elemepts M°'°' = Mf"' = Mf^. The approximation 
of neglecting the interference term was also made by Penn et altB The Coulomb matrix elements for scattering 
between equal and opposite spins are denoted by M°'°' and M'^'^ . We use M^^T = ]\/[li ^nd M'^^ = M^'^ . For the 
calculation, we define two parameters for the average matrix element squared and for the ratio of M^^^ and M^^: 



AfTT 



MTi 



The simplified equations used in the calculations of the scattering rates are given in Eqs. 
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FIG. 1. Illustration of the monochromatic 2PPE process with initial state Ei, intermediate state E2 and final state £3. A 
first photon excites an electron from an initial level Ei in the range between Ep and Ep — hi' into a level £2- The population 
f{E2,t) depends on the temporal pulse shape of the exciting laser and is time-dependent due to electron-electron interaction 
and transport of electrons out of the optically excited region into the bulk. A second photon excites an electron with energy E2 
into a state Es above the vacuum energy Svac, from which it can contribute to the 2PPE intensity via I^^^^{E3,t) oc f{E2,t). 
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FIG. 2. LAPW density of states (DOS) used as input in the calculation of the electron-electron scattering rates in Eqs. 
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FIG. 3. Calculated 2PPE intensity and relaxation time of the excited electron distribution for Cu for photon energy 
hv = 3.3 eV. Curve a shows the 2PPE intensity if no scattering is present and reflects the distribution of optically ex- 
cited electrons. Curve b gives the result if only scattering out of the intermediate level is kept in Eq. (13). Curve c is the result 
if scattering into the intermediate state (secondary electron effect) is also included. Curve d represents the case when the effect 
of transport is also taken into account. The relaxation time when only scattering out of the intermediate state is kept (curve 
6) is a single-electron lifetime and can be compared with the lifetime predicted by Fermi-liquid theory, shown in curve e. The 
other relaxation times are effective relaxation times of the distribution of excited electrons. 
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FIG. 4. Calculated (a) and measured (b) 2PPE intensity and relaxation time for Cu for photon energies hv — 3.0 and 
3.3 eV. Calculations include secondary electrons and transport. Note the dependence of the relaxation time on photon energy, 
especially in the region of the peak in the intensity {E — Ef = 0.7 — 1.3 eV). The minimum in the relaxation time corresponds 
to the peak in the intensity and the maximum in the, ipiaxation time corresponds to the d-band threshold, as indicated by the 
dotted lines, in agreement with several experiments.a"t3 
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FIG. 5. Spin-averaged relaxation time of the distribution and ratio of relaxation times for majority and minority electrons. 
Results labelled owi are single-electron lifetimes. Results {out, in) refer to relaxation times of the distribution including secondary 
electrons. Results {out, in, transport) also include transport. The same average Coulomb matrix element M — 0.8 eV and 
|MfT/Mfi| = 1 is used. 
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FIG. 6. Schematic view of the equal pulse correlation set-up for time resolved 2PPE. 
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FIG. 7. Comparison of experimental relaxation time results for Ag and the three transition metals, Fe, Co, and Ni. 
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FIG. 8. Experimental results for the spin-resolved relaxation time r-f and tj and the ratio r|/r| of majority to minority 
relaxation time. 
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FIG. 9. Experimental and theoretical results for the spin-averaged relaxation time of the distribution. Calculations include 
secondary electron effects. Curve a shows results using Coulomb matrix element M — 0.8 eV and |M^^/M^^| = 1 for the 
various transition metals. Curve 6 gives results for the same M, but |M^^/M^-^| = 0.5. Results using different values of M for 
the various transition metals and |A4^^^/A^^^| — 0.5 are shown in curve c. 
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FIG. 10. Experimental and theoretical results for the ratio r|/r| of the relaxation time of the distribution. The labels are 
as in Fig. H and refer to the same parameters. 
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FIG. 11. Calculated occupation function for different intermediate state energies (normalized to 1 at the maximum) and fit 
by a function describing exponential decay. The calculation is for Cu with hf — 3.3 eV and a pump laser of 70 fs duration. 
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